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Abstract 

>- ■ We determine an approximate expression for the 0(a^) contribution \2 to the kernel of 
the BFKL equation, which includes all collinear and anticollinear singular contributions. 
This is derived using recent results on the relation between the GLAP and BFKL kernels 
(including running-coupling effects to all orders) and on small-x factorization schemes. We 
present the result in various schemes, relevant both for applications to the BFKL equation 
^ ■ and to small- x evolution of parton distributions. 
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1. The BFKL equation beyond next-to-leading order 

Higher order calculations in perturbative QCD, both at fixed order and at the resummed 
level, are playing an increasingly important role in precision collider physics [1]. Fixed- 
order and resummed results pose important constraints on each other. On the one hand, 
fixed-order computations determine an infinite number of coefficients of resummed ex- 
pressions: in the soft (e.g. large-x) limit [2,3], the resummed hard cross-section can be 
determined to any desired logarithmic order by a finite fixed-order computation of suit- 
ably high accuracy, while at high energy (e.g. small x) there is a duality [4] such that 
a finite-order computation of the GLAP kernel determines an all-order resummation of 
the BFKL kernel and conversely. On the other hand, available resummed results de- 
termine partly higher fixed-order expressions and allow an approximate reconstruction of 
their form: in fact, approximate determinations of the next-to-next-to- leading order GLAP 
splitting functions [5] , which played a useful role phenomenological until the exact expres- 
sions became available[6], are significantly constrained by knowledge of the large- ir and 
small- a; behaviour from resummed results. 

The interplay of finite-order and resummed results is especially interesting for the high- 
energy (small- a;) limit of hard cross-sections, which behave as genuine two-scale processes. 
The dependence on the hard scale (henceforth Q 2 , for definiteness) and the energy scale 
(actually, the dimensionless ratio, x for definiteness, that controls the energy dependence) 
are governed by a pair of evolution equations, the GLAP and BFKL equations, respectively, 
whose kernels are related by a duality relation which was recently shown to hold to all 
orders at the running-coupling level [7]. This duality relation determines the resummed 
expansion of either kernel in terms of the fixed-order expansion of the other. Hence, even 
when the symmetry between the two scales is broken by the running of the coupling (and/or 
by kinematics) knowledge of either of the two kernels at fixed orders enables an all-order 
resummation of the other kernel. 

These results have been mainly used to perform small-x resummation of GLAP evolu- 
tion (see refs. [8,9] and refs. therein, and refs. [10,11] for phenomenological applications), 
i.e. to learn about higher-order contributions to the GLAP kernel. However, they can also 
be used to determine higher-order contributions to the BFKL kernel: indeed, they provided 
a nontrivial check on the next-to-leading order determination of the BFKL kernel [12,13]. 
Their use to determine corrections to the BFKL kernel beyond next-to-leading order, which 
are hitherto unknown, has been hampered by two difficulties. First, whereas some next- 
to-next-to leading order duality relations have been worked out some time ago [14], they 
didn't include the full next-to-leading order running of the coupling, nor the informa- 
tion contained in the then unknown NNLO GLAP kernel. In fact, the inclusion of these 
contributions is quite hard if the running-coupling duality relations are solved by brute 
force, even using computer algebra as in Ref. [14]. Second, beyond leading order both 
fixed-order and resummed results, or, equivalently, both the BFKL kernel and the GLAP 
kernel depend on a choice of factorization scheme. Factorization schemes at small- £ were 
determined explicitly in refs. [15,16] and their implications for the GLAP-BFKL duality 
were worked out in refs. [17,10], but only up to next-to-leading order. 

Recent results solve both difficulties. In ref. [7] a general method has been developed 
which allows an efficient determination of duality relations by purely algebraic techniques, 
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with full inclusion of the running of the coupling to any desired order. Thanks to the 
full computation of the GLAP splitting functions to next-to-next-to leading order, [6] a 
computation of all the singular contributions to the next-to-next-to leading order BFKL 
kernel X2 is now possible. Furthemore in ref. [18] small- a; scheme changes have been 
discussed to all orders. Even though in ref. [18] the scheme change required for the BFKL- 
GLAP is determined explicitly only up to NLO, it turns out that its determination up to 
NNLO is possible from available results, as we shall discuss below, at least for the terms 
which affect the singular contributions to the NNLO BFKL kernel. The purpose of the 
present paper is to perform these computations, and use them to construct an approximate 
form of the \2 kernel. 

A determination of the BFKL kernel to NNLO is of theoretical and phenomenological 
interest for various reasons. It is well-known that NLO corrections to the BFKL kernel 
are large and in fact change completely its qualitative shape. The determination of the 
NNLO contribution is thus motivated not only by the slow convergence of the perturbative 
expansion of the BFKL kernel, but also by the expectation that (because of the alternating 
sign of the dominant contributions) the NNLO approximation has a minimum like the LO, 
and thus better stability properties than the NLO. Also, it is unclear whether a direct 
extraction of the NNLO BFKL kernel from the high-energy behaviour of parton-parton 
scattering amplitudes analogous to the NLO computation of Ref. [12] is feasible, because 
it is unclear whether beyond NLO some form of "reggeization" holds, i.e. whether the 
exchange of an effective colorless multigluon state is universal [19,20]. If this is not the 
case, a derivation of the NNLO BFKL kernel from high-energy factorization may be the 
only viable option. 

This paper will be organized as follows. In the next section we will summarize the 
formalism of Ref. [7] for the algebraic resolution of duality relations between the GLAP 
and BFKL kernel, and describe specifically its application to the extraction of all available 
information on the NNLO BFKL kernel from the known NNLO GLAP result. Then in 
Section 3 we will turn to the issue of scheme dependence at small x at NNLO: after 
summarizing the general results of ref. [18], we will describe its application to NNLO 
duality. In Section 4 we will finally determine explicitly the approximate NNLO X2 kernel 
in various relevant factorization schemes, discuss its features and estimate the accuracy 
of our approximation. Technical results on higher-order dualities and on the so-called Qq 
scheme at NNLO are collected in the appendices. 



2. The GLAP-BFKL duality 

Duality is the statement that the solutions to the GLAP and BFKL equations coincide up 
to higher-twist corrections provided their kernels are suitably matched. Its consequence is 
that the leading-twist part of each kernel is determined by the other kernel. This result is 
straightforward to establish at fixed coupling [21,4,22], but rather more subtle when the 
coupling runs [23,7]. Here we summarize the main results, while referring to Ref. [7] for a 
more comprehensive treatment. 

We discuss evolution equations for a parton distribution G{x, Q 2 ), which can be 
thought of simply as the gluon density, or as an eigenvector of a two-by-two evolution 
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matrix in the singlet sector. The kinematic variables x and Q 2 can be thought of as the 
standard DIS variables, or more generally the perturbative scale and a dimensionless scale 
ratio such that < x < 1. We will not consider the dependence on other kinematic vari- 
ables, such as transverse momentum and rapidity, i.e. we will consider standard parton 
distributions, for which ordinary collinear factorization applies. We will therefore limit 
ourselves to angular-averaged quantities at the leading-twist level. We express the parton 
density as a function of the logs of the relevant kinematic variables: 

G = G(£,f); C = log-, t = log%, (2.1) 
and define the Mellin transform with respect to either (or both) of the kinematic variables: 

/•OO 

G(N,t)= / dfe-^Gfof), (2.2) 
Jo 

/oo 
dte- Mt G(^t). (2.3) 
-oo 

Note that, by slight abuse of notation, we denote with the same symbol the parton distri- 
butions G(N, £), G(£, M), and G(N, M), although they are of course different functions of 
the respective arguments. 

The GLAP equation and BFKL equations express respectively the t or £ dependence 
of G. They take the form 

j t G(N,t)=^a s (t),N)G(N,t), (2.4) 
^G(£,M) = *(&., M)G(£,M), (2.5) 

where a s (t) is the running coupling, which upon Mellin transformation becomes the oper- 
ator a 8 , obtained by replacing t — > — in the expression for a s (t). For example, at the 
leading-log level, where a s (t) = a 8 /(l + Poa s t), 

ot s = — — q-, (2-6) 

1 - /3oa S gM 

with a s = a s (0). 

At fixed coupling it is easy to show that if the kernels x an d 7 are related by 

N = x(a s , 1 (a s ,N)), (2.7) 
M = 1 {a s , X {a s ,M)), (2.8) 

then the BFKL and GLAP equation admit the same solution. This relation is straightfor- 
ward to derive [24,22] from the observation that the leading-twist behaviour of G(iV, M) 
is determined by its pole in the (M, N) plane, and that eqs. (2.7) and (2.8) express the 
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position of this pole. This duality maps the expansion of x(ct s , M) in powers of a s at fixed 
M onto the expansion of 7(ct s , A) in powers of a s at fixed a s /N, and the expansion of 
7(a s , A) in powers of a s at fixed A onto the expansion of x(a s ,M) in powers of a s at 
fixed a s /M. 

At the running-coupling level, duality relations can be derived order by order by 
solving eqs. (2.4) and (2.5) perturbatively and matching the respective solutions [14]. That 
this is possible beyond next-to-leading order is highly nontrivial, and it can be proven to 
all orders using operator methods [7]. Namely, one shows that at the running-coupling 
level, the BFKL and GLAP solutions coincide if the respective kernels, viewed as operators 
in (M, A) space, are the inverse of each other when acting on physical states, i.e. such 
that if 

MG(N, M) = 7 (d s , N)G(N, M), (2.9) 

then 

NG(N, M) = X (a s , M)G{N, M), (2.10) 

and conversely. Note that the conditions eq. (2.9) and (2.10) should be viewed as conditions 
on the action of the operators 7(o; s , A) and x(d s ,M) on physical states: specifically, 
they are not just the Mellin transforms of eqs. (2. 4), (2. 5), which depend on a boundary 
condition, but rather, they express a property of the leading-twist Green function of 
perturbative GLAP or BFKL evolution [7]. 

The inversion is nontrivial because the operators involved do not commute. Indeed, 
we can start with eq. (2.9) and construct x(a s , M) as a series in M such that 

X(a s , 7 (d a , N))G(N, M) = NG(N, M). (2.11) 

Because a s and A commute, x(a s ,M) coincides with the inverse function of 7 i.e. the 
fixed-coupling (or naive) dual eq. (2.7): it is the same power series in M. However, because 
7(qi s , A) and M do not commute, we cannot use eq. (2.9) to identify 7(ct s , A) with M in 
eq. (2.11), and thus obtain the desired equation (2.10). 

The problem can be formalized as follows: given an operator equation of the form 



qG = pG, (2.12) 
and given a function f(q), determine the function g such that using eq. (2.12) one gets 



f(q)G = g(p)G. (2.13) 

It is easy to see that when p and q do not commute, the functions / and g do not coincide: 
in fact they can be determined in terms of each other using the expansion of the Baker- 
Campbell-Hausdorff formula [25] to lowest nontrivial order: 

M)G = {fip) - \f"{p) [P, q] + • • • }G, (2.14) 

so that 

9(P) = m-kf"(p)1p,q] + ---- (2-15) 
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Expressions up to third order are derived in ref. [7], and higher-order expressions (up to 
fifth order) are given in Appendix A. 
Hence, if we identify 

p = M, q = >y(a a ,N), (2.16) 

and then choose / = x, we may use eq. (2.15) (and its higher-order generalizations) to 
determine x(a a , 7(d s , N)) on the l.h.s. of eq. (2.11) as a function of M, which we then 
identify with the sought-for BFKL kernel x(d s , M). The computation of x(«s, M) is thus 
reduced by eq. (2.14) to the determination of commutators, which in our case can all be 
obtained from the basic commutator 

[a~\M] = -p + a s poPi + (2.17) 
where the QCD beta function is given by 

P{a s ) = -p a 2 s {l + p 1 a s ) + .... (2.18) 
Using this commutator in eq. (2.14), to lowest nontrivial order we get 



NG(N, M) = 



x (a s , M) + l/3 d^ 7 (« s , N) x"(& s , M) 



G(N,M) + Q(a 2 s ), (2.19) 



where primes denote derivatives with respect to the second argument (M or N) . 

In order for the r.h.s. of eq. (2.19) to provide us with an expression for x(a s , M) we 
must eliminate the residual N dependence in it. This can be done in two steps. First, 
we can back-substitute the lower-order expressions for N given by eq. (2.19) in its higher- 
order terms: for example, eq. (2.19) tells us that we can just replace the leading order 
expression x(d s , M) for the N dependence of the next-to- leading order 0(/3o) term, up to 
0{a 2 s ) corrections: 



NG(N, M) = 



x(& 8 , M) + I/3 a2_^ 7 ( ds? x(a s , M)) x"(a s , M) 



G(N, M) + 0(a 2 a ). 



(2.20) 

Of course, beyond next-to-leading order this back-substitution becomes nontrivial, and it 
must be performed by keeping into account the commutation properties of the operators 
involved. 

The dependence of the result eq. (2.20) on 7(a s , x(d s , M)) and its derivatives can be 
finally expressed using the duality relation eq. (2.8) in terms of M and of derivatives of the 
fixed-coupling dual x °f 7 : e -g- differentiating eq. (2.8) with respect to a s or with respect 
to M gives respectively 

- 7 (a s ,x(« s ,M)) + 7 , (a s ,x(« s ,M)) T — x{a s ,M) = 0, (2.21) 



oa s oa s 
7 / (a s ,x(« s ,M))x'(a s ,M) = l, (2.22) 

and so on. Using this, the right-hand side of eq. (2.20) gives us finally an expression for 
x(ct!s, M) in terms of x(a s , M), i.e. for the running-coupling dual BFKL kernel in terms of 
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the naive dual one, which in turn can be determined from the anomalous dimension 7 using 
fixed-coupling duality eq. (2.7). To lowest nontrivial order this gives the well-known [26,24] 
result 

X (a s , M) = x(a s , M) - \(3 a s x{a s , M) *^ 5 '^ + 0(a 2 s ), (2.23) 

where we made use of the fact that the kernels start at order d s , so <9x(a s , M)/da s = 
a~ 1 x(&s, M) + 0(a s ). Equation (2.23) expresses the BFKL kernel x which is dual to the 
given GLAP anomalous dimension 7, in terms of the fixed-coupling dual x determined 
from 7 using eq. (2.7). 

Analogously, one can construct the GLAP anomalous dimension 7 which is dual to a 
given x m terms of the fixed-coupling dual 7 determined from x using eq. (2.8): 



7 (d s , N) = 7(« s , AO - bo&sXi&s, l(& 8 , V) ?}?''^*'"^ + 0{a% (2.24) 

2 [x(as,l(as,N))] 



where x(d s , 7(d s , A)) can be further expressed in terms of 7(& S ,A) and its derivatives 
differentiating the duality relation between x and 7, analogously to eqs. (2. 21), (2. 22). It is 
important to note that 7 and x are n °t the fixed-coupling dual of each other: rather each 
of them is the fixed-coupling dual of (respectively) x an d 7, which are related by running- 
coupling duality. It is sometimes convenient [8] (and it will be useful for our discussion 
of scheme changes) to view x as an effective x kernel: namely, once 7 is determined from 
a given x using eq. (2.23), we also define x which is the fixed-coupling dual of this 7 
eq. (2.23). 

In ref. [7] we have used this approach to derive the running-coupling corrections to the 
small- ;r resummation of the GLAP kernel which one obtains from BFKL, i.e. essentially 
eq. (2.24), up to next-to-next-to-leading order. In the sequel of this paper, we will use it to 
derive the running coupling contributions to the next-to-next-to leading order BFKL kernel 
obtained starting from the GLAP kernel, i.e. essentially eq. (2.23). Various higher-order 
running-coupling duality relations are derived and collected in Appendix A. 

Because running-coupling corrections to duality are given as a series in a s of terms 
each of which is a function of the fixed-coupling dual expression, also at the running- 
coupling level duality maps the expansion of x(ct s ,M) in powers of a s at fixed M onto 
the expansion of 7(ck s , N) in powers of a s at fixed a s /N, and conversely. This means that 
knowledge of 7 up to next-to-next-to leading order allows us to determine all collinear 
(M ~ 0) singular contributions to x U P to next-to-next-to leading level. Specifically, if we 
expand 

X (a s , M) = a sXo (M) + a 2 sX i(M) + ..., (2.25) 



and then 



»W = ^ + ^ + ..- (2-26) 



for some coefficients Cij, we can determine the first three orders of the expansion eq. (2.26) 
of Xi{M) for all i, i.e. Cij for all i and j = —i — 1, —i, —i + 1. Furthermore, the symmetry 
properties of x allow us to determine its expansion about M = 1 from knowledge of the 
coefficients of the expansion about M = 0. This procedure requires some care in the 
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Figure 1: Comparison of the exact and approximate expressions of the leading-order 
and next-to-leading-order contributions to the BFKL kernels Xo(M) and ^i(M). The 
approximate expressions (given explicitly in appendix D) are obtained by symmetrizing 
the third-order collinear expansion about M = of the kernel, as discussed in section 4. 
The upper two plots are the leading-order xo, while the lower are the next-to-leading- 
order xi' m each case on the left we compare the exact (dot-dash) and approximate 
(solid red) expressions as a function of M, while on the right we show the relative differ- 
ence (exact — approximate) /exact. See section 4 for a full discussion of the construction 
of the approximation and its scheme dependence. 

treatment of the running of the coupling, which affects the way the symmetry is realized, 
and it will be discussed in sect. 4. In the specific case of X2 5 knowledge of 7 up to NNLO 
allows us thus to determine all the singular contributions to X2 at M = and M = 1, i.e. 
all the collinear and "anticollinear" singularities respectively. 

Because xo an d Xi are known exactly, we can test the accuracy of the approximation to 
them which is obtained by retaining only the first three terms of their expansion about M = 
and M = 1 (and subtracting the double-counting between the two expansions). This 
comparison is displayed in figure 1 for the leading and next-to-leading order contributions 
Xo and Xi to the BFKL kernel, where it is seen that the approximation is exceedingly 
good. The percentage accuracy of the approximate expression is better than 0.8% at 
leading order and better than 1.5% at next-to- leading order. 

This motivates us to consider the construction of a similar approximation to X2- Before 
doing this, however, we must turn to issues of factorization scheme dependence. 
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3. Factorization schemes at small x 

Running-coupling duality, reviewed in the previous section, implies that, up to higher- 
twist corrections, the solution to a GLAP-like equation (2.4) can also be obtained as the 
solution to a BFKL-like equation (2.5), and it tells us how the evolution kernels of the two 
equations are related. However, as is well known, beyond leading order, both the GLAP 
and BFKL kernels are only defined up to a choice of factorization scheme. Namely, if the 
normalization of t) is redefined by a subleading function Z(a s , N) = 1 + 0(a s ), the 
evolution kernel beyond leading order changes: hence, the kernel is uniquely defined only 
once the normalization is fixed. 

In general, of course, factorization schemes can mix the singlet quark and gluon with 
each other. Here, however, we will only consider factorization-scheme changes of the single 
parton distribution t) which enters both the GLAP and BFKL equations: because 
only one eigenvector of the anomalous dimension matrix has leading small- a; singularities, 
the BFKL equation is one-dimensional, and mixing is irrelevant for duality. The general 
matching of scheme changes which have effect at small x with those which do involve 
operator mixing (relevant for phenomenology) is discussed in detail in refs. [17,10]. The 
most general factorization-scheme change is then obtained by redefining 

G'(t, N) = Z sx (a 3 , ^) Z lx (a s , N) G(t, N), (3.1) 

where the large- £ and small-x scheme change functions have respectively the form 

Z lx (a a , N) = 1 + a a Z\ x (N) + 0(a 2 s ) (3.2) 
Z sx (a., f ) = Zl (^) + a s Z] x + O(aJ), (3.3) 

with the constraint that Z® x (0) = 1. Upon scheme change, the leading-order contributions 
in the expansion of the BFKL and GLAP kernel in powers of a s remain unchanged, while 
higher-order terms are modified. It is interesting to observe that the lowest-order nontrivial 
small- £ scheme change Z® x (jf) (which affects the NLO BFKL kernel Xi(iV) in eq. (2.25)) 
amounts to a leading-order redefinition of the normalization of gluon, i.e. it starts at order 
a® in the expansion of Z in powers of a s at fixed jf- [17]. 

A further related complication is due to the fact that the BFKL equation is, in its 
most general form, given for a /cj_-dependent parton distribution, and its form eq. (2.5) is 
obtained from angular averaging of this /c^-dependent parton distribution. This angular 
averaging leads to an "unintegrated" parton distribution Q(N,t), which is the derivative 
of the usual parton distribution G(N,t): Q(N,t) = -^G(N,t) i.e., in Mellin space, 

G(N, M) = MG(N, M) . (3.4) 

Because of eq. (2.17), the evolution kernels x(d s , M) and Xi(a s , M) for the BFKL equations 
satisfied by respectively G(M, N) and G(M, N) do not coincide, and are related by 

X (a s , M) = M X i(a s , M)M~ 1 . (3.5) 
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Hence, if the BFKL kernel is determined for the unintegrated gluon distribution Q(N,t), 
while the GLAP anomalous dimension refers to the standard integrated quantity G(N,t), 
there are extra contributions to the duality relation, due to the non-commutativity of M 
and a s in eq. (3.5). The relation between integrated and unintegrated parton distributions 
eq. (3.4) can be viewed as a scheme change by using the (iV, t) space form of the GLAP 
equation (2.4) to relate Q to G, namely: 

g(N,t) = 1 (a s (t),N)G(N,t). (3.6) 

Note that this is not a proper scheme change eq. (3.3) because it does not reduce to the 
identity in the limit a s — > 0. Note also that just as the evolution kernels for Q and G do 
not coincide, similarly the anomalous dimension for their evolution differ: in fact, from 
eq. (3.4) and the GLAP equation (2.4) one finds that Q satisfies 

j t G{N,t) = lu (a s (t),N)g(N,t), (3.7) 

where 

7«M*), AO = 7 (« s (t), N) + j t In [y(a s (t), N)] . (3.8) 

Henceforth, we shall denote by 7 the anomalous dimension at the integrated level, and 
by x the kernel at the unintegrated level, and by 7^ and Xi the unintegrated anomalous 
dimension and integrated kernel. 

Let us now come to the scheme choices which are relevant for the explicit determi- 
nation of j(a s ,N) and x(a s ,M). The normalization of the parton distribution which 
appears in the GLAP equation (2.4) is fixed by the standard factorization of collinear sin- 
gularities [27], and a choice of subtraction prescription such as e.g. dimensional regulariza- 
tion and the MS prescription. This defines anomalous dimensions in the MS factorization 
scheme. Duality then implicitly defines a corresponding factorization scheme for the BFKL 
equation. However, the direct computation of the next-to-leading order BFKL kernel is 
based on the determination [12,13,28] of the gluon Green function in the high-energy limit. 
The extraction of the large energy behavior of the gluon Green function, and the resum- 
mation of its large-energy logs through the BFKL equation are based (explicitly [13] or 
implicitly [12,28]) on a factorization of cross-sections in terms of a high-energy parton dis- 
tribution (the so-called k± factorization [29]) which is compatible with the usual collinear 
factorization, but differs from it by a computable scheme change. 

This is due to the fact that the usual parton distribution which enters the collinear- 
factorized GLAP equation, and the gluon density which enters the fc^-factorization formula 
are normalized differently. This means that, even though the gluon Green-function itself 
is computed in the MS scheme, the evolution kernel extracted from it corresponds to a 
scheme which is not MS, because it describes evolution of a quantity which differs from 
the MS parton distribution by a normalization factor, i.e., it can be obtained from the MS 
parton distribution by a suitable scheme-change function Z sx eq. (3.3). This scheme-change 
function defines the so-called Qo factorization scheme [16]. Furthermore, the quantity 
which naturally enters /c^-factorization formulae is the unintegrated parton distribution 
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Q(N,t), so Qq scheme results are usually given for this quantity, though, of course, Qo 
scheme results can also be given for the integrated parton distribution by using eq. (3.5) 
and the corresponding relation eq. (3.7) between anomalous dimensions. 

The normalization mismatch between k± factorization and collinear factorization, and 
thus the precise definition of the Q scheme, has been determined in ref. [29] at leading 
nontrivial order, i.e. at the level of Z® x eq. (3.3), which affects the definition of xi m 
the expansion eq. (2.25) of the BFKL kernel and therefor its dual GLAP anomalous 
dimension 7 up to next-to- leading order (order 7 SS ) in the expansion of 7(ct s , N) in powers 



of a s at fixed jf: 



7 (a s , N) = 7s (^) + a s7ss (^) + . . . . (3.9) 



The scheme change at the next order (relevant for % 2 and 7 SSS ) has been recently derived 
in ref. [18]. 1 

The main result of ref. [18], which we shall use in what follows, is an expression (proven 
up to NNLO, but conjectured to hold in general) which relates the t dependence of the 
integrated parton distributions G(N, t) (as defined in standard collinear factorization) in 
the MS and Qo scheme. This relation is expressed in terms of the BFKL kernel for the 
unintegrated distribution Q(N,t). Specifically, the t dependence can be written in terms 
of a saddle-point evolution factor E(t,to), a running-coupling duality correction J\f(N,t), 
and a normalization factor 1Z(to) which is characteristic of the way minimal subtraction 
with dimensional regularization is defined by continuation of the anomalous dimensions in 
d dimensions. 

The saddle-point evolution factor is obtained by solving the running-coupling BFKL 
equation for the unintegrated distribution in the saddle-point approximation: this can 
be shown [30,26,23] to lead to evolution driven by the anomalous dimension 7 u (a! s (t), N), 
obtained from the unintegrated BFKL kernel using naive (fixed-coupling) duality eq. (2.7), 
but with a s = a s (t), 



E(t, t ) = exp 



%(a 8 (t'),N)dt' 



to 



(3.10) 



The running-coupling correction to duality discussed in the previous section can be 
combined with the factor eq. (3.6) which relates the integrated and unintegrated parton 
distributions. This gives 



Af(N,t) -y(a a (t ),N) 



exp j jf [7«M*'). N) - 7uM0, N )] dt ' 



AT(N,t ) <y(a s (t),N) 

= exp I jf [ 7 (a s (0, AT) -%(a a (t'),N)] dt' j . 



(3.11) 



1 Note that in ref. [18] (and elsewhere) the scheme change effected by Z° x is referred to as 
a leading-log x scheme change, essentially because it is a leading-order redefinition of the gluon 
normalization, whereas we will call it a NLO scheme change, because it affects the NLO kernel 
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Note that lim a s _> J\f(N, t) = 1 so this can be viewed as a scheme change in the proper 
sense. 

Finally, the normalization factor 1Z is due to the fact that in the MS scheme the 
anomalous dimension is defined as the residue of the simple e pole in the partonic cross 
section, which in turn is given by 



7(<* B ,iV,e) = 1 
(3(a 8 ,e) a s e 

1 



1 PM 



a s e 



+ ...) (^(a 8 ,N)+ej(a 8 ,N) + ...) 



7 (a s ,iV)-^7(a s ,iV) + ...), 



(3.12) 



where in MS /3(a s , e) = a s e + f3(a s ) is the <i-dimensional j3 function, 7(a s , e) is the anoma- 
lous dimension obtained using duality from a <i-dimensional BFKL kernel, which has been 
expanded as 

7 (a s , N, e) = 7 (a a , N) + £ 7 (a s , N) + 0(e 2 ). (3.13) 

It follows that the anomalous dimension determined from duality differs from the MS 
result through the terms beyond the first on the r.h.s. of eq. (3.12), and thus MS evolution 
requires an additional factor 



K(N, t) 
K(N,t ) 



— exp 



/3 a(t')i(a(t'),N)dt' + 0(a 2 a ) 



t 



(3.14) 



where we have used (3(a s ) = —floa 2 . + . . .. 

Combining all these factors, the result of ref. [18] takes the form: 



G Q « (AT, t) = Af(N, t)E(t, t )K(N, t )G MS (N, t ) 



(3.15) 



This equation gives the scale dependence of the parton distribution in either scheme in 
terms of a boundary condition determined in the other scheme: therefore, it fully spec- 
ifies both the relation between the two schemes, and the scale dependence in either of 
them. Letting t = to in eq. (3.15) immediately gives the relation between G®°(N,t) and 
G MS (N,t), through the function 



R{N,t) =J\f(N,t)H(N,t). 



(3.16) 



The scale dependence of the parton distribution in the Qo scheme is found by keeping 
to fixed in eq. (3.15) and varying t: thus 



= exp 



<y(a s (t'),N)dt' 



to 



G Qo (N,t ). 



(3.17) 
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X I Ai 

Figure 2: Schematic relation between BFKL kernels \ an d GLAP anomalous dimensions 
7 in various schemes. Horizontal lines denote duality relations while vertical lines denote 
scheme transformations. The diagonal lines express the identities eq. (3.21) and (3.22). 



So the integrated parton distribution in the Qo scheme evolves with the anomalous dimen- 
sion 7(iV, t) which is related to the starting unintegrated BFKL kernel by running-coupling 
duality combined with the transformation to the integrated level. 
The scale dependence in the MS scheme is instead given by 



G MS (N,t) 



n(N,t ) 

K(N,t) 
K{N,t ) 
K(N,t) 



E(t,t )G MS (N,t 



exp 



%(a s (t'),N)dt' 



(3.18) 



G MS (iV,t ), 



namely, the integrated parton distribution in the MS scheme evolves with an anomalous 
dimension which is closely related to the fixed-coupling dual j u (N,t) of the starting (un- 
integrated) BFKL kernel, and only differs from it through the scale dependence of the 1Z 
factor eq. (3.14). In this respect the MS scheme is remarkably simple. 

In fact it is possible, and useful, to define an auxiliary scheme, MS , which differs 
from MS by a factor 1Z(N, t): namely 

G m * (N, t) = K(N, t)G m (N, t) 



(3.19) 



In this MS scheme the parton distribution evolves with the naive dual anomalous dimen- 
sion: 



G MS (N,h 



exp 



%(a s (t),N)dt 



to 



G MS (N,t ). 



(3.20) 
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The relations between different quantities in various schemes, which will be computed 
in the next section, are summarized in figure 2. In the figure, horizontal lines denote 
duality: either at the running-coupling level, relating x to 7, or at the fixed-coupling, 
relating 7 to \. Vertical lines denote relations between schemes, specifically the Qo, MS 
and MS schemes. Equation (3.20) means that the anomalous dimension in the MS scheme 
coincides with the naive dual of the Qo scheme BFKL kernel (at the unintegrated level): 

^°(a s ,N) = 1 m *(a s ,N), (3.21) 

and thus, by duality, also that 

xF\a s ,M) = X Q0 (a s ,M), (3.22) 

where x^ s (a s ,M) is the naive dual of the standard GLAP anomalous dimension in 
the MS scheme, while x < ^°(q; s ,M) is the kernel for the BFKL equation satisfied by the 
unintegrated parton distribution in the Qo scheme. This further implies that if we interpret 
X^ s (a s , M) as an operator by letting a s — > a s , and order it canonically (i.e. in the same 
way as x^°(d s , M)), then it is related by running-coupling duality to 7 Qo (a s , N). These 
relations are denoted by diagonal lines in the figure, and will turn out to be useful in 
relating results obtained in various schemes in the next section. 



4. The NNLO BFKL kernel 

The construction of the NNLO BFKL kernel proceeds in three steps. First, we collect the 
results discussed in the previous sections to get expressions of the BFKL kernels which are 
related by either running-coupling duality (x(ct s , M)) or fixed-coupling duality (x(a s , M)) 
to a given GLAP kernel 7(0^, N) to NNLO in various schemes. Then, we use the known 
expression of the GLAP anomalous dimensions up to next-to-next-to leading order to 
determine the first three coefficients of the expansion eq. (2.26) of the leading Xo(&s, M), 
next-to-leading xi(&s,M), and next-to-next-to-leading x 2 (a: s ,M) order BFKL kernels in 
powers of M about M = 0. Finally, we use the underlying symmetry of the high-energy 
gluon emission diagrams to determine the corresponding coefficients of the expansion of 
Xo(a s , M), Xi(d s , M) and X2(« s , M), about M = 1. 

Because all duality relations are expressed in terms of naive (fixed-coupling) dual 
kernels, our starting point is the determination of the BFKL kernel x(a s ,M) which is 
obtained from the NNLO GLAP anomalous dimension using the fixed-coupling duality 
equation (2.8). Expanding x(et 8 , M) in powers of a s at fixed a s /M as 

X(a s , M) = x s (a s /M) + a s x ss (a s /M) + ... (4.1) 

we have 

7o(x.(a./M)) = — , (4.2) 
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Xss(a s M) = 7T — — — — , (4.3) 



Xsss 

(a s /M) 



1 



%( Xs (a s /M)) 
1 



l2 ( Xs (a s /M)) + 1 ' l (x s (a s /M)) Xss (a s /M) 

(4.4) 



+ ^{Xs(a s /M))(x ss (a s /M)) 2 

Expanding out both 7^ and x s i, i = 1,2, 3 up to third order in their respective ar- 
guments, eqs. (4.2)-(4.4) determine the coefficients of these expansions in terms of each 
other. Starting with the expansion of 7 

7o (iV) = ^ + 9o,o + 9o,iN + 0(iV 2 ), (4.5) 

li(N) = ^+gi,o + 0(N), (4.6) 

72(A0 = 9J + + 0(N°), (4.7) 



we get 



Xo(M) = ^ + 0(M 2 ), (4.8) 
X\(M) = + ^ + |^ + O(M), (4.9) 



XUM) 



(9o,-i) 9o,i + 9o-i(9o,o) , 9i,o9o,-i + 9o,o9i,-i , 92 -1 r>(M°\ (4 10) 
M 3 ~ M 2 ' ' M 



where the superscript i reminds us that all these quantities are defined at the integrated 
level. The values of the coefficients gij up to NNLO in the MS scheme can be extracted 
from the known GLAP anomalous dimensions [31], and are listed in appendix B. Using 
these expressions of g^, eqs. (4.8)-(4.10) give the first three terms in the expansion of the 

naive dual kernel x m the same scheme: this corresponds to relation between 7 MS and 
X MS in figure 2. Note that the vanishing of leading singularities of 7 at NLO and NNLO, 
(7i _2 = (72,-3 = 0, implies the well-known vanishing of the constant and linear term in the 
LO BFKL kernel (4.8). 

Starting from the naive dual kernel x\ e Q- (4.8)- (4. 10) we can determine the BFKL 
kernel in the same approximation in various factorization schemes. For the sake of com- 
parison with direct diagrammatic computations, the Qo scheme is most relevant, because a 
direct computation using minimal subtraction at the level of the BFKL equation gives the 
kernel in this scheme: indeed, available expressions [12,13,28] of the next-to-leading kernel 
Xi provide the result in the Qo scheme. Using duality, this kernel can be obtained from 
the MS GLAP anomalous dimension by exploiting eq. (3.21) (see figure 2), which states 

that x^° is the naive dual of the (integrated) anomalous dimension 7 MS . Hence, we can 
obtain by first transforming the anomalous dimension to the MS scheme through the 
1Z scheme change eq. (3.14) and then using naive duality eqs. (4.8)-(4.10). 
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The scheme change which takes the anomalous dimension to MS can be determined, 
using eq. (3.14), from the expression of the GLAP splitting functions in d = 4 — 2s di- 
mensions 7(AT, e), as explained in Appendix C. The MS anomalous dimension is given 
by 

7 MS* {N) = 9±-l + Q + ^ + 0{N 2^ (4 n) 

tF^^+Slo + OW, (4.12) 
lF S \N) = 9 -^ + ^ + 0(N% (4.13) 



where all coefficients are the same as in the MS eqs. (4.5)-(4.7) except 

9i,o = 0i,o - Po9o,o , 92,-i = 92-i -Po9i,-i - Pl'go-i. (4.14) 

The coefficients <7o,o and (70,-1 are determined by the <i-dimensional LO splitting func- 
tions; their values are given in Appendix C. The coefficient <7i 5 _i is determined by the 
(i-dimensional NLO GLAP kernel which is as yet unknown. This is the only term con- 
tributing to the singular part of the NNLO BFKL kernel that we have not been able to 
calculate explicitly: we will estimate below the uncertainty related to our ignorance of this 
contribution to the scheme change. 

The BFKL kernel in the Qo scheme is thus given in terms of the naive dual kernel 
eq. (4.8)-(4.10) by 

X$"(M) = x (M), (4.15) 
X? (M) = x\(M) + O(M), (4.16) 
v Qo / M \ _ (ffo,-i) 2 ffo,i + 9o ,-i(go,o) 2 , <7l,QgQ,-l +90, 09l, -1 92,-1 n(M°VA 1 7\ 



We can also determine the BFKL kernel in the MS scheme, by observing that it is 
related by running-coupling duality to the MS anomalous dimension; note that this now 
gives the kernel at the integrated level (figure 2 again). Hence, we add to eqs. (4.8)-(4.10) 
the running-coupling corrections discussed in Sect. 2, and explicitly given in appendix A, 
eqs. (A.13)-(A.15) in terms of the naive dual expressions computed above. Defining 

A m X (M) = X P(M) - xUM), (4.18) 
and expanding about M = we get 

A m X o(M) = (4.19) 
A M5 Xi(M)=/3 ^ (4.20) 

A^S-v (M\ - a ( 3 90,090,-l ggi -1 2 ff2,-2 \ , na 2 90,-1 , a a 90,-1 (A r>, x 
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Using eqs.(4.19)-(4.21) in eq. (4.18) we get the BFKL kernel in the MS scheme at 
the integrated level. The expression at the unintegrated level, to be compared to the 
unintegrated Qo scheme result of eqs. (4.15)-(4.17) can be obtained using using eq. (3.4), 
which implies 

±Q{t M) = (a sX o + Ajxi + & 3 sX 2 - [&., M] g - [a 2 s , M] g + 0{a A 8 )) M), 

(4.22) 

so the kernels x\ a t the integrated level and Xi at the unintegrated level are related by 

Xo = X > 

Xi=Xi- Po^, (4.23) 

X2=X2-/3oA§-2/3 §. 

So far, we have exploited the information from the GLAP anomalous dimension to 
determine the first few terms in the expansion of the BFKL kernel about M = (the 
collinear region). However, we can further use the underlying symmetry of the BFKL 
kernel to determine the corresponding terms in the expansion of the BFKL kernel about 
M = 1 (the anticollinear region). At LO and NLO this leads to an approximation to the 
BFKL kernel which, when tested on the known LO and NLO kernels, turns out to be 
accurate to better than 2% in the whole physical region < M < 1, as shown in figure 1. 
Indeed, the underlying Feynman diagrams which determine the unintegrated BFKL kernel 
are symmetric upon the exchange of the incoming and outgoing gluon. This means that the 
dimensionless BFKL kernel K(a s , k 2 , Q 2 ), related to x(d s , M) by Mellin transformation 

x(* s , M ) = l°°^ {%) M K K k2 > Q 2 ) > ( 4 - 24 ) 

is symmetric upon the interchange of the virtualities of the incoming gluons Q 2 and k 2 : 

-^K(a s , k 2 , Q 2 ) = ^K(a s , Q 2 , k 2 ). (4.25) 

This symmetry, in turn, implies that the BFKL kernel x(a s ,M) is symmetric upon the 
interchange M <-»■ 1 - M [12]. 

However, the kernel is symmetric only if one chooses a symmetric argument for the 
running coupling, and if the iV-Mellin eq. (2.2) does not break the symmetry between 
the two scales Q 2 and \i 2 which enter the definition eq. (2.3) of the M-Mellin transform. 
In deep-inelastic scattering both symmetries are broken: the argument of the running 
coupling is Q 2 , and £ = ln(s/Q 2 ). Hence, the BFKL kernel which is obtained from the 
GLAP anomalous dimensions incorporates these symmetry-breaking effects. The symme- 
try breaking must be undone before the symmetry can be exploited. After symmetrizing 
one can revert to nonsymmetric variables and argument of the coupling in order to get an 
expression of the dual x which is accurate for all M. 
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A symmetric choice of variables is for example £ = ln(s/ a/ Q 2 k 2 ). As is well 
known [12], the kernel x s (ct s ,M) which corresponds to this choice is related to the kernel 
X DIS (a s , M) which corresponds to DIS variables by the implicit equation 

X s (d s , M) = X DlS (& a , M + \ X s {a s , M)) . (4.26) 

A symmetric choice of argument for the running coupling is such that eq. (4.25) also 
holds when the argument /j 2 of the running coupling a s (fJ. 2 ) is expressed as a function 
of Q 2 and k 2 , y 2 = y 2 (Q 2 ,k 2 ), and y 2 (Q 2 ,k 2 ) = y 2 (k 2 ,Q 2 ). Examples of symmetric 
choices for the running of the coupling are /j 2 = \Q 2 — k 2 \ or /j 2 = Max(Q 2 , k 2 ). Upon 
Mellin transformation, different choices of arguments of the running coupling correspond 
to different orderings of the running-coupling operator. For instance, it is easy to check 
that the Mellin transform of 

K (a s ,k 2 ,Q 2 ) = ^a s (Q 2 ) ^ Q f_ (4.27) 



is 6t s xo{M), where 



X o(M) = -— [,/>(M) + if>(l - M) - 2^(1)] (4.28) 

7T 



is the standard leading-order BFKL kernel, while the Mellin transform of 

K (a s ,k 2 ,Q 2 ) = V as(fc2) |Q2 Q2 fc2 | (4.29) 

is Xo(M)a s . 

Note finally that the integration which takes us from the unintegrated to the integrated 
distribution also breaks the symmetry, as eq. (3.5) explicitly shows: it follows that the 
symmetry which holds at the unintegrated level is broken at the integrated level. 

Because the diagrammatic computation of the BFKL kernel yields the result in the 
Qo scheme at the unintegrated level, once symmetric variables and symmetric running of 
the coupling have been chosen the kernel is symmetric when determined in this scheme. 
Whether the symmetry is preserved or not in other schemes depends of course on the par- 
ticular scheme change. Specifically, it is easy to see that a scheme change through 7Z(N, t) 
preserves the symmetry at leading nontrivial order (i.e. at the order of xi) but not be- 
yond, essentially because it is entirely determined by xo an d its (/-dimensional continuation, 
which is symmetric (see appendix C), whereas a scheme change through Af(N, t) breaks the 
symmetry already at leading nontrivial order, because the first nontrivial running-coupling 
duality correction eq. (2.23) is manifestly not symmetric. 

As discussed already, use of naive duality on the MS GLAP anomalous dimension 
eqs. (4.11)-(4.13) gives us the expansion eqs. (4.15)-(4.17) of the Qo scheme BFKL kernel 
for the unintegrated parton distribution. This result is clearly in DIS variables and can be 
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turned into symmetric variables by expanding out eq. (4.26). Substituting the expansion 
eq. (2.25) of the kernel we get 

X s (d s ,M) = a sXo (M + \a sX s (M) + \a 2 sX \(M)) 

+ a 2 sX i {M + \a sX s (M)) + a^ 2 (M) + 0(aj) (4.30) 
= a sXo (M + \a aX s Q {M)) + a 3 sX 'o {M + \a sX s (M)) \ X \(M) 

+ a 2 sXl (M) + afoi MjxSM + a 3 sX2 (M) + 0(aj). (4.31) 

where on the right-hand side we have dropped the DIS index of eq. (4.26). 

The first term on the right-hand side of eq. (4.31) must be computed by carefully 
keeping operator ordering into account. This can be done by using a technique akin to 
that of ref. [7], summarized in section 2: namely, by computing 



M,|a sX( s ,(M) 



d 



2 



-+■ 



Xo(M+id,x5(M)) = e" — n " 'J'*' Xo(M + A)| A=0 

= Xo (M) + \a sX s {M) X ' G {M) + \ [M, a sX s (M)} X q(M) + \a 2 sX s2 {M) X ' \M) + 0(a s s ) 

= Xo (M) + \a sX s (M) X ' G (M) - \foa 2 sX l(M) X ' \M) + \a 2 sX s2 (M) X '^M) + 0(a 3 s ). 

(4.32) 

Substituting this result in eq. (4.31) and collecting terms of the same order we get 

XS(M) = X ?°(M), (4.33) 
X \{M) = X ?°(M) + |xo Q °(M) X J°(M), (4.34) 
X|(M) = X ?°(M) + ixi Q °(M)x?°(M) + i X o Qo (M)x?°(M) 



2A1 V-^/AO J 1 2- 

+ ko Q °( M ) 2 Xo Q °(^) + U /Qo (M)x Q ° 2 (M) - |/3 x[, /Qh (M)x? (M), (4.35) 



where on the right-hand side we have denotes the result in DIS variables by x& , since in 
practice we will use the Q -scheme result in DIS variables of eqs. (4.15)-(4.17). 

In order to determine the constant term of xi and the simple pole of X2 in symmetric 
variables we have substituted the expansion of xo up to 0(M 2 ). In principle the linear 
term of xi is needed too, but its dependence cancels out in the expression for X2- Using 
the expressions eq. (4.15)-(4.17) for the unintegrated Qo-scheme kernels x^° we finally get 

xi(M) = ?2j=l + 0(M 2 ) , (4.36) 
XKM) = + + ^ + ^ 4- 9o 2 _ lC (3) + O(M) , (4.37) 

y s (M) = ^°-^ 3 _ 3 9oA9o,-i) 2 + Poido-i) 2 (9o,o) 2 9o,-i + ffo,i(ffo,-i) 2 ~ 9o-i9i-i 
V - J 2M 5 2M 4 ' " M 3 

-502,-2+00,001,-1+00,-101,0 , 92,-1 -2/3 (ffo_i) 2 C(3) n , A/f( u , 4 „m 
+ ^ + m +0{M). (4.38) 
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We can now exploit the symmetry of the kernel with symmetric coupling, which implies 
that x(ct s , M) must admit an expansion of the form 



X (a s , M) = a sX s Q (M) + 6l 2 s x\{M) + d* X £(M) + 0(a 3 s ) 

= X S o ym (&s, M) + x? m (« s , M) + X s 2 ym (&s, M) + 0(a 3 s ) 

where x| ym (d s ,M) are the symmetrized functions 

1 1 



(4.39) 
(4.40) 



X s vm (a s ,M)=co,- 1 



ot* 1 ot s 

M 1 — M s 



+ d s c ,o + cq,i [a s M + (1 - M)a 6 



+ c ,2 [a s M 2 + (1 - Mfa s ] + 0(M 3 ), 



(4.41) 



xf m (d s ,M)= Cl ,_ 3 



1 



Q' 



+ 



M 3 (1-M) 3 s _ 



+ Ci _ 2 



+ 



+ Cl -1 

+ 0(M 2 ), 



a °M + 1-M as 



5 = 1,5 



1 



+ 



(1 - M)J 



M 2 (1-M) 2 s _ 
+ a s c 1)0 + ci,i [c^M + (1 - M)a 2 s ] 

(4.42) 

+ O(M ). (4.43) 



a 



Given the expression of Xi Vm e Q- (4.40) it is straightforward to determine the sym- 
metrized kernel when a s is "canonically" ordered to the left, which corresponds to the 
choice of argument of the running coupling a s = a s (Q 2 ): 



a s x s ym (M) = Xo Jm (&s,M), 



(4.44) 



alx? m (M) = X f m (« s , M) - 6^o ^ + Po& 2 s (c ,i + 2c , 2 ) + O(M), (4.45) 



a 



My 

co,-i 



Cl,-1 ,~?,n C l-2 C l-3 



2 ^°(1-M) 2 



Aa 6 s f3 Q 



(1-M) 3 



(1-M) 4 



+ 0(M°). (4.46) 



In x^ ym the symmetry is broken by the running of the coupling only. 

We can finally determine all coefficients Cij in eqs. (4.41)-(4.43) by expanding the 
symmetrized kernel eqs. (4.44)- (4.46) in Laurent series about M = and equating to the 
expansion of the unsymmetrized kernels xf e Q s - (4.36)-(4.38), which is accurate to the 
stated power of M. Because the anticollinear terms with poles at M = 1 in eqs. (4.41)- 
(4.43) are regular in M = 0, the symmetrized x sym have the same M = poles as 
their unsymmetrized counterparts xti an d their coefficients can be read off eq. (4.36)- 
(4.38). However, the anticollinear terms do contribute to all regular contributions in the 
expansion of xf about M = 0. This is why higher-order regular terms must be included 
in the right-hand side of eqs. (4. 41), (4. 42): specifically, symmetric terms up to 0(M 2 ) 
must be included in Xo(M) in order for its expansion to coincide with that of Xo Vm (M) 
up to and including O(M), and terms up to and including O(M) in Xi(M) in order for its 
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0.2 0.4 0.6 0.8 



M 

Figure 3: The approximate expression of the NNLO contribution to the BFKL kernel for 
the unintegrated distribution eq. (4.46) in the Qo scheme, symmetric variables, a s (Q 2 ). 
See appendix D for the values of all coefficients. The uncertainty band is obtained 
varying the unknown scheme — fixing coefficient —5 < gi,~i < 5 (top to bottom). 




M 



Figure 4: The relative uncertainty in the NNLO contribution to the BFKL kernel for 
the unintegrated distribution shown in figure 3, due to the uncertainty in the unknown 
scheme — fixing coefficient gi,-i, here taking the values —5, —1, f , 5 (top to bottom). 
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Figure 5: The full BFKL kernel at leading, next-to-leading, and next-to-next-to leading 
order obtained combining tyhe known expressions for the LO and NLO contributions 
(as shown in figure 2) and our approximate expression (with <?i,-i = 0) for NNLO (as 
shown in figure 3), with a s — > 0.2. All expressions are in the Qo factorization scheme, 
at the unintegrated level with symmetric variables, and a s = a s {Q 2 ). The symmetry 
about M = h is only broken by the argument of the running coupling. 

expansion to coincide with that of x^ m (M) up to O(M ). No addition is necessary for 
X2 because the known coefficients in its expansion about M = are all singular. 

Summarizing, all singular coefficients in eqs. (4.41)-(4.43) can be read off eq. (4.36)- 
(4.38), while for the nonsingular ones we get 



co.o 

Cl,0 



■§0o,- 



92,- 



C0,1 = 0, 
2, 



9o,- 



f + Q7o,-i) 2 C(3) 



§0o,- 



C0,2 : 

ci,i = l(go 



1) 

\2 



<7o,-i#o,o — 9i,-i- 



(4.47) 



Using these results in eqs. (4.44)-(4.46) we get our approximate expression for the BFKL 
kernel up to NNLO, at the unintegrated level in the Qo scheme in symmetric variables, 
with a s = a s (Q 2 ). The LO kernel of course does not depend on either scheme, the choice 
of variables, or the running of the coupling. The NLO kernel corresponds to the widely 
used form of the kernel as given in ref. [12], eq. (14) of that reference. Indeed, it can 
be straightforwardly checked that the Laurent expansion of eq. (4.45) coincides with the 
result of ref. [12] up to and including O(M ). The NNLO kernel eq. (4.46) is a new result. 
The expressions for the kernel in D1S variables can be obtained straightforwardly from 
eqs. (4.44)-(4.46) by inverting eqs. (4.34)-(4.35), and the kernel at the integrated level is 
found using eq. (4.23). The MS scheme expressions are found using eqs. (4.20)-(4.21) in 
eq. (4.18). 
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Figure 6: The BFKL intercept c(a s ) at leading, next-to-leading, and next-to-next-to- 
leading order, obtained from the x as shown in figure 5, with a symmetric argument for 
the running coupling, and evaluated at M = |. 



The full expression for the kernel up to NNLO in the Qo scheme is given in Ap- 
pendix D, both in symmetric and DIS variables, at the unintegrated and integrated level. 
The approximate expressions of the LO and NLO kernels compared to the exact kernels 
in figure 1 correspond to the expressions eq. (4.44)-(4.45) (Qo scheme, unintegrated, sym- 
metric variables, a s (Q 2 )). The approximate expression of X2 is displayed in figure 3, in 
the Qo scheme, in symmetric variables, at the unintegrated level and with the canonical 
argument of the coupling a s (Q 2 ). As discussed above, eq. (4.14) (see also Appendix C) 
one of the coefficients which determine the next-to-next-to-leading order scheme change 
between the MS and MS (and thus Qo), namely <7i,-i, is unknown. This coefficient affects 
the simple (sub-subleading) poles, and has therefore a moderate impact. Noting that all 
the scheme-change coefficients (see Appendix C, eq. (C.7)) are of order one or smaller, 
and indeed all coefficients gij (see Appendix B, eq. (B.2)) are at most of order of a few, 
we estimate the uncertainty related to this coefficient by varying —5 < gi,-i < 5. The 
corresponding uncertainty is displayed in figure 3, and the relative uncertainty in figure 4: 
it is seen to be similar to the uncertainty of a few percent that we expect (on the basis of 
the LO and NLO results of figure 1) to affect our approximate form of \i- 

In figure 5 we display the full NNLO BFKL kernel x(« s , M) = a s xo(M) +a 2 s xi(M) + 
a^X2(M), with the same scheme and variable choices using the exact expressions up to 
NLO and the approximate expression to NNLO, with a s — > 0.2. The slow convergence 
properties of the expansion of the BFKL kernel, driven by the increasingly dominant 
collinear and anti-collinear singularities at M = and M = 1, are very apparent in this 
figure. Although the X2 contribution restores the minimum near M = §, the convergence 
of the expansion in the vicinity of the minimum is still rather slow. 
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This is made more explicit in figure 6, where we plot the intercept c(a s ) as a function 
of a s . This is calculated by using the same expression for \ as in figure 5, but with the 
coupling chosen to be symmetric, so that all the curves are symmetric about M = |, and 
then defining c(a s ) = x(a s , -|), ie the value of \ the stationary point. While the fixed 
order perturbation series is clearly good for very small ct s , say a s ~ 0.05, with a s ~ 0.1 
there are signs that the series has yet to converge. For yet larger values of a s (as would be 
appropriate for phenomenological studies) the results from the fixed order series are clearly 
not very useful, and a resummation of collinear and anti-collinear singularities along the 
lines discussed in Ref. [33,9,22,8] becomes necessary. 

5. Outlook 

In this paper we have presented an approximate determination of the NNLO contri- 
bution to the BFKL kernel. In the process, we have provided a full treatment to this order 
of various issues which affect the determination of the BFKL kernel: the relation between 
the MS and Qo factorization schemes, the duality relations which connect the BFKL kernel 
to the GLAP anomalous dimension, specifically in the presence of running coupling, the 
choice of kinematic variables in the definition of the BFKL kernel, the relation between 
the form of the BFKL kernel and the argument of the running coupling, and the relation 
between BFKL kernels for integrated and unintegrated parton distributions. All these 
issues become rather nontrivial to next-to-next-to leading order, and require full control 
of factorization scheme and running coupling. 

Because the perturbative expansion of the BFKL kernel in both the collinear and 
anticollinear regions is alternating in sign, a knowledge of NNLO corrections is necessary 
for an accurate assessment of the uncertainty involved in a fixed-order determination of the 
kernel: indeed, whereas the qualitative features of the NLO kernel are completely different 
from those of the LO, the NNLO result is qualitatively similar, though we have shown 
that it is quantitatively not so reliable because of the slow convergence of the perturbative 
expansion, even in the central region away from the singularities. 

Fixed-order BFKL kernels have been widely used recently in studies of nonlinear 
(saturation) corrections to the BFKL equation and their phenomenological implications 
for RHIC and the LHC (see e.g. [34] and ref. therein). Also, they are the foundation 
of numerical approaches to the BFKL equation (see [35] and ref. therein), which in turn 
are relevant for Monte Carlo simulations (see e.g. ref. [36] and ref. therein). Because 
of the slow convergence of the perturbative expansion, the determination of the NNLO 
BFKL kernel presented here is useful in assessing the reliability of these calculations. 
Finally, the approximate form of the LO and NLO kernels given here are extremely accurate 
while having a very simple analytic form, and are thus amenable to simple numerical and 
phenomenological implementations. 
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various discussions during the course of this work. S. M. is supported through a SUPA 
graduate studentship. S. F. was partly supported by a PPARC visiting fellowship and a 
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Appendix A. Higher— order duality 

Higher-order duality relations can be obtained by pursuing to higher orders the ex- 
pansion of the Baker-Campbell-Hausdorff equation for a pair of non- commuting operators 
p, q which act in the same way on a state G, according to eq. (2.12). To fifth order we get 



f(q)G(N,M) 



m-im \p,q] + hm [qamh 

+ \f{P)'" [p, [p, q}] + y{p) IV [q, [q, [q,p]}} + 
+ \f{P) IV [q, IP, IP, q}}} ~ \f{P) IV IP, [p, [p, q}}} + 
+ y{p) IV \PA] 2 - lf(P) V [PA] [q, [q,P]} + 



(A.l) 

f ( H\ y I i-i ri\ \ri I ri ■to I I _l_ 

24- 

- hf{p) v [<?, [q, p]] [p, q] - hfi.P) v % q] % % q}\ + 



y{p) V [P, PAW [P,q] ~ hf{p) VI \pA? + O(aj) 



G(N, M) . 



We can now obtain higher-order generalizations of the BFKL-like equation obtained 
starting from a GLAP equation by a suitable identification of p and q. Because, as well- 
known (see e.g. [7]) at the fixed-coupling level the dual of the expansion of x(ct s ,M) in 
powers of a s at fixed M is the expansion of 7(a s , N) in powers of a s at fixed a s /N, it is 
convenient to expand 

7 (d s , AdJ 1 ) = -fsiNa' 1 ) + d^AdJ 1 ) + . . . (A.2) 

We then perform the identification eq. (2.16) of p and q and use eq. (A.l) with f(p) = 
x(d s ,p), where 

X(a, 7 (d, Ad" 1 )) = Ad" 1 . (A.3) 
Namely, x(d s , M) = d7 1 x(d s , M), where x is the naive (fixed-coupling) dual of 7, so that 

X(d, M) = xo(M) + axi(M) + d 2 x 2 (M) + ... . (A.4) 

Using eq. (A.l) we then get 

Na~ l G{N,M) = \x(a,M) - ±[M, 7 ]x"(d, M) - \[M, [M, 7 ]]x'"(d, M)+ 

l (A.5) 
+ -[7, [7, M]]x"'(d, M) + §[M, 7] W, M)}G(A, M) , 

where primes denote derivatives with respect to M, and we define 7 = 7(d s , Na^ 1 ). 
All commutators can now be determined explicitly using the expression 

1 = i - i3 °Wm + ft (- aA a¥ - ^ )2 Mp) + 0(a -» (A ' 6) 
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of the running coupling at the operator level: 



[7, M] = - (NP + Nfo&sPi ) d ^Z N * sl) + 

oNa s 

+ A^^£^ + 0(a»), (A.7) 

:.3> 



[ 7 ,[ 7 ,M]] = 0(^) (A.8) 

d(Nal 



[M,[M, 7 ]] = (iV/3o) 2 ^|^^ + 0(^), (A.9) 



([7,M]) 2 = Wo f ( ^y^ + Q^)- (A.10) 

Substituting the commutators (A.7)-(A.10) in eq. (A. 5) and back-substituting order 
by order the low-order expansion of the enduing equation in the higher-order terms one 
may remove all the N dependence from the right-hand side of eq. (A. 5), with the result: 

Na s G = < x( a s, M) - ^a s (3 ^— — + 



2 



+ a 2 s 



i*2-, A (Qs,7) 1r a x(a s ,M) X "(a s ,M) 



X'(a«»7) 2 x'(« s ,7) 

+ ^ %5?' M S ( i2 ^"^- M )) 3+ (A - n) 

24 X '(a s ,M) V 

- 7x'(« s , M)x"(« s , M) X '"(a s , M) + 3(x'(« s , M)) W,, M)) + 

Identifying the term in curly brackets in eq. (A. 11) with the BFKL kernel 

X(&s, M) = a sXo (M) + a sX i(M) + a 2 s x 2 {M) + ... . (A.12) 

and expanding X as i n e Q- (A. 4), eq (A. 11) gives an order-by-order expression of the 
running-coupling dual in terms of the naive dual. Up to NNLO we get 

Xo = Xo (A.13) 

Xi = Xi-§A)^ (A.14) 
Xo 

Inn XOXO , 
X2 = X2 - 2^0^!—^- + 

Xo 



+ ^° 2 w ( 12 ( *° )3 " u * w> + 3 ( *° )2 x ~ IV °) + 

',11 -T, ?,H?,I 



In AUA1 n A1A0 , In AUAO 

Xo Xo (xo) 

l a 2 Xo fo-r/ ~/" 
(Xo) 



+ |/3o 2 7^ 2x Xo "(Xo) 2 , (A.15) 
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where all x% an d X% are functions of M and the prime denotes differentiation with respect 
to M. 

By inverting this relation (i.e. expressing order by order x i n terms of x), this result 
can be used to determine the running-coupling corrections to the anomalous dimension 7 
determined from a given BFKL kernel x- These corrections up to NNLO were first derived 
in ref. [14], and then reproduced more recently in refs. [18], [7]. However in both cases only 
the leading-order kernel xo and leading-order running of the coupling were included. For 
completeness, we derive them here consistently including the running of the coupling up 
to NLO. 

Namely, we start from a given BFKL kernel x(a s , N) and determine first the GLAP 
anomalous dimension j(a 8 , NaJ 1 ) which is obtained from it using naive duality, which 
we then expand according to eq. (A. 2). We wish to determine the anomalous dimension 
7(a s , Na^ 1 ) which is related through running-coupling duality to the starting x- We 
observe that the naive-duality relations between x an d 7 and between x an d 7 imply that 

N/a s = x(a a , 7(a s , N/ot a )) = x(a s , 7(a s , N/a a )). (A.16) 

Note that in eq. (A.16) it is immaterial whether a s is considered to be an operator or not, 
because [a s ,N] = anyway. Defining 

A X (a s ,M)=x(a s ,M)-x(a s ,M) 

(A. 17) 

A7(d s , Na s 1 ) = 7(d s , Na s 1 ) - 7(d s , Na s ), 
and expanding eq. (A.16) in powers of a s we get 

Axi = -* iA7 " (A.18) 

A% 2 = -x A7 sss - XiA7 ss ' - x'i^Iss + |xo A 7ss + Xo A lss'A-f SSl 

where Axi(M) and A^ii^Na^ 1 ) are respectively the coefficients of the expansion of 
Ax(ct s ,M) and A7(d s , Na^ 1 ) in powers of a s , and on the left-hand side all Axi are 
evaluated as functions Axi(7(ct s , Na^ 1 )). 

Furthermore, inverting eq. (A. 15) we can determine the coefficients Axi of the expan- 
sion of Ax(a s , M) in powers of a s in terms of x (all Xi functions of M): 

* _ in X0X0 
^Xi — 2^0 — — 
Xo 

A _ 1 o a X0X0 1 1 o X0X1 1 o X1X0 i Q X0X0X1 in , 

AX2 - 2P0P1— — + 2^0—^— + Po— ; 2^0 2 (A.19) 

Xo Xo Xo (Xo) 

+ ^o{^ (6 (XoO 3 " lOXoXo'Xo" + 3 (x ) 2 Xo' v ) + \dl 7^ Xo'- 

Equating the right-hand side of each of eqs. (A.18) to the corresponding equa- 
tions (A.19) we determine finally 

A 7ss = -i/3 ^ (A.20) 



Xo 



M=7 S 
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and 

A 7sss = A 7 W + A 7 « + A A 7ss , (A.21) 

where we have defined 

A 7 ^S = jg^ (15 (Xo) 3 " ISXoXoXo + 3 (x ) 2 Xo^) (A.22) 

and 

Ar rill = (§A)XoXo7m + §A)XoXi7" + §A)XoXi7" + §A>XoXo7"7m) „ . • ( A - 23 ) 

The full NNLO running coupling correction eq. (A.21) is given here for the first time. 

In particular, the term A 7 iss, which depends only on xo an d Po, was already derived 
in ref. [14] (eq. (14) of that reference) and later confirmed in ref. [18] [eq. (4.16) of this 
reference, which gives Af(N,t) eq. (3.11)]. However, the terms A 7 sss, which are due to the 
next-to-leading order kernel %i, and the last term on the right-hand side of eq. (A. 20), due 
to the next-to-leading order running of the coupling were never computed before. 



Appendix B. Expansion of the GLAP anomalous dimensions 

The coefficients of the expansion eqs. (4.5)-(4.7) of the leading, next-to-leading and 
next-to-next-to-leading GLAP anomalous dimensions in powers of N are easily determined 
by recalling that 7 is the large eigenvalue of the 2x2 anomalous dimension matrix, given 
by 

1= \ Igg + 111 + V in 99 ~ 7 99 ) 2 + ^Igglgg , (B.l) 



and using the expressions of 7 ^ given in Refs. [31,32,6]. In the MS scheme we get 
9o,-i = Q f 



9o,o - — i^r + {--^ + 3^Cl ) n f 

n _ Cat* , 67C A HCpn; ( C% , C F \ 2 
( IZC F 23C a \ „ 



^ - 1 V 36tt^ ) J 

_ 2C(3)Cl 1643Cj llCj / 4 3C A , C F 547C F , 4 L ; 
91,0 — ^ I" 216tt 2 36 "r" ^54^" "r 18 216tt 2 ^ 4tt 2 Ca / / 

, / 13C F _ 13C| \ 2 
^108tt 2 C a 547r 2 C 2 y 'V 

#2,-2 
#2,-1 

301Cj , 8C(3)C F CU , 35C f Ca _ 28853C F C A _ 2C F C(3) , HC 2 F \ 
81tt 3 3ti- 3 108tt 25927I- 3 37r 3 127T 3 y '*/ 

, / 59C A _ 65C F \ „2 
I 648^ 324T3"J 



~ 2T67F 2 36 ^ ^ 54T 2 ^ 18 2167T 2 ^ 4tt 2 Ca J J 

, 13C F 13C F \ 2 (B - 2) 

^io8tt 2 Ca 547r 2 c 2 ; "7 

C(3)Cj , llCi _ 395Cj / Ci _ 71Cj _ C F C A , 71C F C A \ 
27T 3 72tt 1087.-3 ^36tt 1087T 3 IStt 54tt 3 y "/ 

143C(3)C^ _ 297rCj _ 389Cj 73091(7% / 11C(3)C A _ C\ 
24^ 720 432tt 25927i"3 + ^ 12^ 9tt 

, 301C 2 , , 8C(3)C F C A , 35C p Ca 28853C p Ca 2C p C(3) , 11C 
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Appendix C. The 1Z scheme change to next-to-next-to- leading order. 



The 1Z scheme change [15,37,18] is related to the fact that the MS anomalous dimen- 
sion contains interference terms between collinear poles and the j3 function in d dimen- 
sions. The factorization of collinear singularities for a <i-dimensional partonic cross section 
& which depends on a single dimensionful variable Q 2 can be written as 



o 



(C.l) 



where a s (fx 2 ) is the dimensionless renormalized coupling, «q is the bare coupling, 
<j(°\Q 2 , ao, N, e) is the regularized cross section and cr(^2 , OL s {p?), N, e) is free of collinear 
singularities. The factorization scale is fx 2 and 7(a s ,iV, e) and (3(a s ,e) are respectively 
the (i-dimensional anomalous dimension and beta function. The former is defined as the 
Mellin transform of the d-dimensional Altarelli-Parisi splitting function. The latter is 
given by 

P(a 8 ,e) = a s e + (3(a s ) (C.2) 
in terms of the usual four-dimensional (3 function 



-p a 2 8 (l + a s (3i + ...). 



(C.3) 



The MS anomalous dimension is the residue of the simple pole in e in the integrand 
of the exponential in eq. (C.l), namely 



7 MS (« s ,iV) = Res £ 



a s 7( a s, N, e) 



= l(a s ,N) 



P{a s ,e) 
(3(a s ) 



a* 



a. 



(C.4) 



where the various coefficients are defined through the Taylor expansion 

7 (a s , N, e) = 7 (a s , N) + £ 7 (a s , N) + e 2 ^y(N) + .... 



(C.5) 



The MS anomalous dimension thus receives two different classes of contributions: pure 
collinear singularities and interference terms between the ^-dependent anomalous dimen- 
sion and the poles arising from the expansion of the (i-dimensional (3- function. In particular 
up to next-to-next-to leading order we have 



7o - 7o 



7 MS 
MS 



7i + A)7o 



(C.6) 



72*~ = 12 + A)/3i7o + /3 7o + Mi- 

The 1Z scheme change, which appears in the relation between MS and Qo schemes 
eq. (3.15), takes us from the MS scheme to the MS scheme where the anomalous dimension 
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is simply given by 7(a s , N, 0), i.e. 70 (AT), 71 (AT), etc. Thus in order to compute 7 MS we 
have to subtract from the MS anomalous dimension eq. (C.6) the contributions coming from 
the interference between the <i-dimensional kernel and the (3 function, which are in turn 
determined from the knowledge of the <i-dimensional anomalous dimension "f(a a , N,e). 

The d dimensional leading order splitting functions have been known for a long time, 
at least for x < 1 [38]: 



1 



P qq (x,e) = C F — 
P qg (x,e) = C F 



(1-xY 
1 



1 + x' 



P gq (x, e) — Tj 



R 



(l-x) e 
1 

(1 -x) s 



1 — X 
l + (l-x) 2 



e(l-x) 



+ a qq (x,e)8(l - x) , 



— ex 



x 



1 - 2x- 



1 



x 



(C.7) 



P gg (x,e) = 2C A - ( —^ e 



1-e 

x 1 — x 



1-x 



+ 



X 



+ x(l — x) 



+ a gg (x,e)8(l - x) . 



The end-point contribution a qq (a gg ) can be extracted from any process with collinear 
radiation from incoming quarks, such as Drell-Yan, or gluons, such as Higgs production 
from gluon fusion: the 0(a s ) coefficient of the 5(1 — x) provides a determination of the 
end-point term in the splitting function after factoring a simple e pole and the Born cross 
section (and a factor of two when there are two incoming partons) . Using the known NLO 
corrections for Drell-Yan [39] and Higgs [40] production we get 



a qq{ £ ) = Cf 



+ - +e 4- 



7T 



+ 0(£ 2 ), 



/ \ 2C^ 11C A - 4n f T R 2 2 
a gg( e ) = — + g-^ en 2 +0(e z ). 



(C.8) 



The simple e pole cancels against that coming from the expansion of (1 — x)~( 1+£ ^ = 
-6(1 — x) + . . . in the splitting functions P qq and P gg: thereby providing a check of the 
result. 

Using the splitting functions (C.7)-(C.8) we can now determine the coefficients of the 
expansion in powers of iV of the large eigenvalue of the ^-dependent GLAP anomalous 
dimension matrix: we find 



7o(A0 = 9 ~^ + 9 -^ 1 +9o,o + 0(N), 
70 (AO = ^ + O(N ) 



(C.9) 



with 



90,-2 = 0, go-i = 0, 



7T 



67 7 rif 
^ = _ T27"8TV 9o '- 1 = ~12 



(CIO) 
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The next-to- leading order d dimensional splitting function is not available, though in 
principle it could be extracted from <i-dimensional splitting amplitudes [41]. However, the 
1Z scheme change has been determined long ago [15] to the level of the next-to-leading 
N = singularities, i.e. for 7 ss (a s /iV) eq. (3.9). This corresponds to the 0(e) correction 
to the leading N = singularities, because [18] 



7f S = 7* 



7 ^ S = 7ss + /3o7 s , 



(C.ll) 



where we have used a notation similar to that of eq. (C.5) but for the expansion eq. (3.9) 
of the anomalous dimension. Hence, if we let 

71 = ^ + ^ + ^ + 0(N°) . (C.12) 

the coefficients <7i,-3 and #1,-2 can be extracted using eq. (C.ll) from the scheme change 
ofref. [15]. 

Equation (C.ll) implies that 7 S is left unaffected by the scheme change, so it follows 
from eq. (C.6) that </o,-i = #0,-1 = 9i,-3 = thereby confirming the result of eq. (C.10). 

Also, the scheme change [15] of 7 SS starts at O (a s (jf) 3 ^j 
Collecting everything we get 

<7i,_ 3 = 0; <7i,_ 2 = 0, (C.14) 

while the sub-sub leading coefficient (71,-1 remains undetermined: it would require knowl- 
edge of the O(e) correction to the simple iV-pole contribution to 71 (a s , N, e). 

Summarizing, the MS* anomalous dimension is given in terms of the coefficients 
eq. (B.2) of the expansion of the MS anomalous dimension and of the scheme change 
coefficients eqs. (C. 10), (C.14) by 

7<r* = 9j ^ + goo + go,iN + o(N*), 

7P* = ^+ 2i,o + 0(N), (C.15) 

where 

9i,o = 9i,o - Pogo,o , ^ 
92,-1 = 92,-1 - (3q9i,-i ~ 0q9q,-i • 



30 



Appendix D. The BFKL kernel in the Q scheme 

In this appendix we give explicit expressions for our approximation to the NNLO 
BFKL kernel in the Qo factorization scheme. The kernel for evolution of the unintegrated 
distribution, with the argument of the strong coupling chosen as a s (Q 2 ), and the symmetric 
choice of kinematic variables eq. (4.24) is given by 

Xo(M) = ^(±+ (T^y - 1 - M(l - M)); (D.l) 

^t M \__C\l 1 , 1 \ i C'a ( 11C A n f C F n f \( l l n 

XiyiVl)— ^ [jfw -T (i_ M )3 ) I" n \ 127r 6tt ^ 3ttC a /IFT (l-M) 2 J 

_ CUfl 1 I /13CVn/ _ 23C A n,M fj_ , 1_\ 

7T ^0 (1-M) 2 "r I 18tt 2 36tt 2 Ailf T 1-Mj 

/ CjC(3) C(3)Cj , llCj 395Cj Ca»/ 71C A n/ C F ry /pj 2 ) 

" r I 7T 2 2tt 2 "t" 72 IO871- 2 " r 36 108tt 2 18 V • V 

71C F n/ ^ Ca\ _ Ca / _ 11C A _ nf_ , Cf«j \ 

" r 54tt 2 "r" ^0 7T >/ "r" 2tt 2 ti" I 12tt 6tt 3ttCa I 

X2(M) = (ivp- + (i_M)s) - 2^ ( - - ^ + + A)) (a?* + (1-M)0 

4_ 4« c i i _l ca 17 _ nc A _ hj_ _i_ cw ^ 2 

~T ^0 7T 2 (1-M)4 7T LI 127T 67!" 37!-Ca / 

I CU f C J- n / I CFn J HC F n f C A ir | 67Ca ^ ( l3C F n f 23C A "/ \1 

7T I 9C 3 7T 18C 2 7T 36CA7T 6 3671- / I 187T 2 367T 2 J J 

( _L I 1 \ I oCja (a I IICa i n / 2CVn^ \ 1 

' \M 3 "•" (1-M) 3 J -r z tt PO^PO + 67r "I" 37r 37rC - A ; (i_M) 3 

, [ _ 1 / gjC(3) , llgj _ 395Cj Cjn/ _ 71C A n f _ C A C F n f 71C A C F n f \ 

+ I 2 I 2?T 3 + 72-K 1087T 3 + 367T 1087T 3 18tT ~+~ 54tt 3 J 

4_ ( — HCa n f I Pf"/ ~\ (13Cpnf_ _ 23C A n f \ 

\ 127T 67T 371-Ca / V 187T 2 367T 2 J 

Ca / _ 2C(3)Cj 16430% _ llCj 43C A n f C F n f _ 547C F n / 13C F n^ 
"T" 7T I tF^ r 216tt 2 36 Mtt^ r 18 216tt 2 <~ 1087t 2 Ca 

C|n / _ 13C|n| A,(JZ_ 4- 4- ffl^ll ( 1 4- 1 ^ 

4tt 2 C*a 547r 2 C 2 ~^ P°\12tt 81tv ) ) 2^ J KM* (1-M) 2 J 

r _ 143C(3)C% _ 297rC% _ 389Cj 7309 1C| _ llC%C(3)n/ _ C A n f 
L 24tt 3 720 4327T 2592tt 3 12tt 3 9tt 

301Cjnj 8C(3)CACFn/ SSCaCf^/ . 59C A n^ _ 28853CACFn/ 

8l7T 3 37T 3 1087T 6487T 3 25927T 3 

3lr3 324^ + —V2^ P0^1-l + Po 12 - l ^ tt 2 J Im + T^mJ' 

(D.3) 

Substituting the numerical values of the Casimirs Ca = 3 and Cf = f we get 

Xo (M) = f ( JL + _i_ _ 1 _ M (l - M)), (D.4) 

V1 fMl - A i 1 - 4- J^^-L- - _ 1 103n/ /J_ _j 

A1K 1V1 ) - 2tt 2 VM 3 ^ (1-M) 3 ^ \ 47T 2 187T 2 / M 2 I 2tt 2 9tt 2 ) (1-M) 2 1087T 2 Im f 1 — M/ 

I 11 I nf 143 , 47n / 27C(3) , M (_ 33 , n/ \ 
^ 8 "r" 108 T21F 7 ~T 162tt 2 "r 2tt 2 "t - V 4T 7 "r 2^^ ' 

(D.5) 
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v n (M\ - JJ-(-L_ -4- 1 A I ( 99 I n f\ 1 I ( 495 5nM i 
X2^W J — 27r 3 l,M5 ^ (l-M)B J Utt 3 ^ ttSJ M 4 ^ Utt 3 tt 3 ; (1-M) 4 

i r 1167 i 35n/ n / _ _9_1 1 , r 1893 , 23n/ _ 7n^ _ jjj l 

LT^f 3 " T^f 3 " 1087I- 3 27rJ M 3 " ' L16t^ ~9tz t ' 36t7^ 2ttJ (1-M) 3 

, [ 1653 , 37771/ _ 5n^ , _99_ , W _ 243C(3) l l 

" r Ll6?r 3 "t" 432ti- 3 648tt 3 " r 16tt ^ 24tt 4tt 3 J M 2 

, r 1653 5049 , 88171/ 933n/ _ 211 "/ (r^ c\ 

L 16tt 3 8(33-2n/)7r 3 144tt 3 8(33-2n/)7r 3 648tt 3 

_ 19r^ i _99_ i W _ 243C(3) l 1 

4(33-2n/)7r 3 "r 16tt ^ 24tt 4tt 3 J (1-M)' 2 

, [ 121 _ llw/ | n j | 73091 _ 6125n/ llw? _ 389 _ llgi-i 
L 192 144 432 ^ 96tt 3 648tt 3 1944tt 3 16tt 4tt 

8n/ gi-m/ _ 87tt _ 1683C(3) , 457n/C(3) -| / j_ 1 \ 

27tt ~t 6vr 80 8^ 1 108tt 3 JVM l-MJ' 

The expression of the NNLO kernel in DIS variables can be obtained by inverting 
eqs. (4.30)-(4.31), with the result 

DIS _ ( 9_ _i_ 1167 _ lire/ , n / \ 1 , ( 863 , 33 _ 54C(3) , 235n/ 

A2 ^ 2tt "t 16tt^ 12~k^ 108tt 3 / M 3 " \T6tt^ 4tt tt 3 h 432tt 3 

2n/ _ , /121 I 73091 _ 87tt _ 389 _ Hgi-i _ 1287C(3) 

9tt 6487I" 3 / M 2 ' V 192 96tt 3 80 16tt 4tt 8tt 3 

_ lln/ 8n/ 6125n/ ffi-m/ 133n/C(3) nf_ llnf x l 

144 27ti- 6487T 3 6tt 108tt 3 432 1944-zr 3 / M 

i _ 54 I 1683 _ 31n/ \ 1 

(1-M) 5 7T 3 ^"slF 3 " "l^W (1-M) 4 

_i_ / _9_ _i_ 1893 , 65n/ _ 7n 2 x j_ 

"r" V 2tt "t" 16tt 3 "r" 127T 3 36W (1-M) 3 

, / 33 , 2137 _ 135C(3) 7n f 381 In f _ 211n 2 f x j_ 

V8?r 16tt 3 2tt 3 ^ 36?r ^ 432tt 3 648ti- 3 ) [1-M) 2 

4- ('121 i 73091 _ 389 _ 87tt _ Hgi-i _ 1287C(3) _ lln/ 8n/ 
^ 192 96tt 3 16tt 80 4vr 87F 3 144 27ti- 

_ 6125n f gi-m/ 133n/C(3) lln/ nj x 1 

648tt 3 6ti- "r" 108tt 3 19447T 3 432/ 1 — M* 



(D.7) 



Finally, the kernel for the evolution of the integrated parton density can be obtained 
from the unintegrated one through eq. (4.23). Using DIS kinematics the difference at 
NNLO is given by 

YUM) - Y U (M) - ( - ^!LL + ^L)^-+ 3(-33+2n/) i 

X2\ m ) X2\ m ) — \ 18-7T 3 277r 3 / M 3 27T 3 (1-M) 3 

12393-4938n/+206n 2 . l -15147+1182n/-16n 2 1 
^ 648tt 3 M 2 "I 108tt 3 (1-M) 2 

+ 37^3 [ - 703890 + 37974n f + 284n 2 f + 29403/t 2 (D.8) 
- 158471/7T 2 - 12nj7r 2 + 96228C(3) - 5832n / C(3)] ^ 



M 

324tt 3 1-M ' 



24543+705n / - 109^ 
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